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We study the de Haas- van Alphen (dHvA) oscillation in two-dimensional and quasi-two-dimensional 
systems. We give a general formula of the dHvA oscillation in two-dimensional multi-band sys- 
tems. By using this formula, the dHvA oscillation and its temperature-dependence for the two- 
band system are shown. By introducing the interlayer hopping t z , we examine the crossover from 
the two-dimension, where the oscillation of the chemical potential plays an important role in the 
magnetization oscillation, to the three-dimension, where the oscillation of the chemical potential 
can be neglected as is well know as the Lifshitz and Kosevich formula. The crossover is seen at 
4t z ~ 8tabH I 'cj>o, where a and b are lattice constants, 4>o is the flux quantum and 8t is the width 
of the total energy band. We also study the dHvA oscillation in quasi-two-dimensional magnetic 
breakdown systems. The quantum interference oscillations such as /3-a oscillation as well as the 
fundamental oscillations are suppressed by the interlayer hopping t z , while the /3+a oscillation 
gradually increases as t z increases and it has a maximum at t z /t ~ 0.025. This interesting depen- 
dence on the dimensionality can be observed in the quasi-two-dimensional organic conductors with 
uniaxial pressure. 



PACS numbers: 71.18.+y, 71.70.-d, 71.70.Di 
I. INTRODUCTION 

It is well known that the magnetization (M) period- 
ically oscillates as a function of the inverse magnetic 
field (1/H) with frequencies of the extremes of the cross- 
section area of the Fermi surface. This oscillation, known 
as the de Haas- van Alphen (dHvA) oscillation, is used to 
study the shape of the Fermi surface of the conductors 
and the experiments have been fitted by using the Lif- 
shitz and Kosevich (LK) formulaau. From the fitting 
of the temperature-dependence of the amplitude of the 
dHvA oscillation, we can obtain the cyclotron effective 
mass (m) for each cyclotron orbits. 

Recently, a revival of interest in dHvA oscillations in 
two-dimensional and quasi-two-dimensional conductors 
arises, since it has been noticed that thepLK formula 
can not be applied in the two-dimensional3~II3 and quasi- 
two-dimensionaltj conductors. Because of the Landau 
quantization of the two-dimensional kinetic energy in the 
magnetic field, both electron number (N) and chemical 
potential (/Lt) cannot be fixed simultaneously as the mag- 
netic field is changed. In the actual situation, N is fixed, 
and /i oscillates periodically as a function of 1/H. In 
the three-dimensional jcanductors the oscillation of /i is 
shown to be very smalltln. Therefore the LK formula ob- 
tained under the condition of fixed /x (grand canonical 
ensemble with treating fi independent of H) is justified 
to be used in the system of the fixed N (canonical en- 
semble). However, in the two-dimensional system, the 
oscillation of /x as a function of 1/H cannot be neglected. 
The saw-tooth waveform of the dHvA oscillation in the 



fixed N case is inverted from that in the fixed /x case in 
the single-band casefl. 

It has been shown that there appear the quantum 
interference oscillations with the frequencies such as 
fp ~ fa and fp + f a in the magnetization when N 
is fixed in magnjetitpbijeakdown systems of the tight- 
binding electronaffi0iiJiL3, in the, magnetic-breakdown 
systems treated— semi-classicallyflO'Ej, and in multi- 
band systemstj EiMlii The anomalous spin-splitting 
effects par e— . also shown to exist in the fixed N 
systemscIrEHl. The quantum interference oscillations have 
been confirmed by the experiments of the dHvA ef- 
fect in various two-dimensional and quasi-two dimen- 
sional systems *uph as the quasi-two-dimensional or- 
ganic conductorsEl (/c-(BED^TF) 2 Cu(NCS) 2 B'P 1 ajcid 
a-(BEDT-TTF) 2 KHg(SCN\Jp), Sr 2 RuO Ja~0, 

GaAs/AlAs heterointepfaceEj and In^Gai-^As single 
quantum well structureO. 

From. .numerical calculations for the single-band 
systemsEKil, it has been shown that the temperature- 
dependence of the amplitude of the dHvA oscillation can 
not be described by the LK formula. The conventional 
fitting of the LK formula may not give the cyclotron ef- 
fective mass in the two-dimension. 

In this paper we give a general formula for the dHvA 
oscillation in two-dimensional multi-band systems. From 
this formula, the oscillatory parts of the chemical poten- 
tial and the magnetization are calculated numerically in 
the two-band system. Moreover, by increasing the inter- 
layer hopping t Z} we study how the dHvA oscillation in 
the two-dimensional system change to those in the three- 



1 



dimensional system. 

In the next section we give the general formula for the 
dHvA oscillation. In section 111 we show the difference 



of dHvA oscillations in cases of the fixed fj, and the fixed 
TV in the two-dimensional multi-band system. In Section 
[V the crossover from two-dimension to three-dimension 
is studied. In section [v| the quantum interference os- 
cillation are studied in two-dimensional and quasi-two- 
dimensional magnetic breakdown systems. 



II. FORMULATION 

The thermodynamic potential il is calculated as a func- 
tion of the chemical potential fx, magnetic field H and 
temperature T as 



n(fi,H) = -Tlog 



^exp 

. N.E 



E-pN" 
f . 



= -T / deN {e, H) log 1 + exp - 



T 



(1) 



where E is the total energy of the system, N is the num- 
ber of electrons and No(e,H) is the density of states in 
the presence of the magnetic field H . The magnetization 
M is obtained in the grand canonical ensemble as 



dQ(n, H) 
dH 



(2) 



As mentioned by Lifshitz and Kosevichlj], /i is treated 
as the independent variable of H in the above equation. 
This is not the case in the actual situation. 

The electron number N is calculated in the grand 
canonical ensemble as 



N(p.H) 



dn((i,H) 



(3) 



If we take fi to be independent of H, N oscillates as a 
function of H . In order to study the system with constant 
N, i.e. canonical ensemble, we should consider Eq.(||) as 
the equation giving the chemical potential p(N, H) as a 
function of N, H and T, 



N 



dSl(ti(N,H),H) 
dn(N,H) ■ 



(4) 



In this case fj, oscillates as a function of H . Then we get 
the Hclmholtz free energy as 



F(N, H) = Sl(n(N, H), H) + fj,(N, H)N 



(5) 



The magnetization should be calculated in the canonical 
ensemble as 



M N (N,H) = - 



dF(N, H) 
dH ' 



(6) 



Using Eq.(|) we get 



M N (N, H) 



dn(fi{N,H),H) _ dn(N,H) 

dH dH 
( dn(j*,H) 
\ dH 



M^(N,H),H) 



N 



»=H(N,H) 



(7) 



We see that the magnetization calculated in the canon- 
ical ensemble has the same form as the magnetization 
calculated in the grand canonical ensemble, but we have 
to take into account the oscillation of /i in the system of 
the fixed electron number. This is the general formula 
independent of the dimension of the system. 



III. TWO-DIMENSIONAL ELECTRONS 

A. multi-band system 

In this section we study the two-dimensional multi- 
band electrons in the magnetic field. This model can be 
applied to the multi-band materials such as Sr2Ru04 and 
In x Gai_ x As. The density of states in the magnetic field 
is given by the sum of delta functions, 



j n=0cr=±l/2 



(8) 



where j is the band suffix, pj is the density of states for 
each spin in the j band at H = 0, 



eH 



m j 



(9) 



is the cyclotron frequency with the cyclotron mass rrij, 

£jna = Uj(n + -) + gj<rfj, B H + Aj, (10) 

is the Landau level, Aj is the energy of the band bottom 
in zero magnetic field, gj is the electron ^-factor for the 
j band and [Ib — eh/2moCo is the Bohr magnetron with 
electron mass mo and velocity of light cq- In this section, 
we set h = cq — ks = 1. The electron number at H = 
is given by 



N 



E E / ^ e = 2 E ft (a*o - a ;)> ( n ) 



j CT=±l/2 



where /xo is the chemical potential at H = 0. 

Applying the Poisson foripulaa in Eq. ([!]), thermody- 
namic potential is written asE3 



n(p,H) = n (fi,H) + ^H), 



(12) 



where 
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Cl ((i,H)=-T / deJ2 T 



Pi 



3 <r=±l/2 

xlog(l + ex P r-^-^ g - e n (13) 



and 



= Y2 E^ 



E E x x 
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where 



Rr.jr 



2ir 2 rT/LUj 



sinh(27r 2 rT/ujjY 
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and 



figjr = COS r7r 



m 7 

9 j = 9 3 — ■ 

m 
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(15) 
(16) 
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In the above we have used the identity 



sin a; 7rcosM/_B) 
r dx = 1 



1 + e 



BsinhtTr/S) 
(-l)" +1 e- An 



(18) 



where A = (// - Aj - gjap B H)/T, B = Uj/(2icrT), and 
we have neglected the summation term over n, since we 
are interested in low temperature properties (A 3> 1). 
At low temperature (T <C fi — Aj), we get 

n (n,H)n~J2 E P3^-^~9j^^bH) 2 . (19) 



3 ff=±l/2 



We obtain 



dH ~ A 



(20) 



From this equation, it is found that d£lo(p, H)/dH does 
not make the oscillation of the magnetization as a func- 
tion of 1/77. Thus, in grand canonical ensemble, (l(p, H) 
gives the oscillatory part of the magnetization M^fj,, H) 
as 



M^(p,H) = - 



dH 



RTjrRs.jr 



7j , !' 



x sin 2nr 



H 2 



(21) 



where 



27re/j = 2nrrij(fi — Aj), 



(22) 



is the area of the Fermi surface at H = in the j band. 
Conventionally, Rr,jr is known as the temperature re- 
duction factor of the LK formula and is used to estimate 
the cyclotron effective mass (rrij) from the T-dependence 
of the amplitude of the oscillation. The electron number 
is obtained from Eq. (||) as 



N(fx,H) 



dfl dfl 
dp dp 



+ Y / —RT,3rR S . 

r=l 



3 r 



x sin 2~nr 



M ~ A? ■ 1 



(23) 



Here we consider the fixed TV case. As we wrote in the 
previous section, the above equation should be taken as 
the equation for p(N,H). With the use of Eq. (|lT|) , we 
get, 



fi(N,H) 



where 
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(24) 



(25) 



and p(N, H) = p(N, H) — po is the oscillatory part of the 
chemical potential. This is a transcendental equation for 
p(N,H) and it is difficult to give p(N,H) by solving 
Eq. (|24|) analytically except for simple cases such as the 
single-band case at T = that we will discuss below. 
After we obtain p(N,H), we can calculate Mn(N,H) 
( Eq . (Q)). From the second equation of Eq. (0) and Eq. 



M N {N, H) 



dQ Q (p,H) dn(p,H) 



dH 



dH 



(26) 



It is found from Eq. ( |20| ) that (3fi (/i, H)/dH is indepen- 
dent of /i. Thus, the first term of of the right side of Eq. 
( p6| ) does not contribute to the oscillation of Mn(N, H) 
as a function oi\/H. As a result, we get 

= M fl (p(N,H),H), (27) 
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where M N (N, H) is the oscillatory part p£ M N (N, H). 

Recently, Alexandrov and BratkovskjO asserted that 
the oscillations of the magnetization also come from flo 
through the oscillation of p. Their result for the free 
energy (Eq. (13) in their paperE3) is formally correct 
but they did not take account of the magnetic field de- 
pendence of the f a , which cannot be neglected in two- 
dimensional systems. As a result their analysis of the 
dHvA oscillation for fixed N system (canonical ensemble) 
is insufficient and their conclusions on Fourier transform 
intensities are incorrect. 

When all of the band-bottom energies are the same 
(Aj = A), we can derive the simple relation between 
p(N,H) and M(N,H). Since £l(p,H) rapidly oscillates 
as a function of (p — A) /H in this case and the envelope 
depends on H slowly, we can approximate 



dCl(n,H) dCl((j,,H)/j,(N,H)-A 



dH 



Using Eq. m) we get 



dp 



H 



goQyg) =N | an 

dp dp 



and using Eqs. (11) and ( |l9| ) we get 
dQ(fi,H) 



-2pp{N, H). 



(28) 



(29) 



(30) 



Then Eq. (p7|) is rewritten as 



M N (N, H) 



Pt ioKN,H) 1 



H 



if Aj = A for all j, 



A | ji(N,H) 
Mo Mo 



(31) 



where \p(N,H)/p \ <C 1 is used. Thus the oscillatory 
part of the magnetization is proportional to the oscil- 
latory part of the chemical potential, if all the band- 
bottom energies are same. The similar relation between 
Mn {N, H) and jl(N,H) in the single-band system with 
electron (or hole) reservoirs has been, obtained indepe, 
dently by Mineev and ChampelEja, Itskovskyei al 
and Grigorievtil. However, when the band-bottom en- 
ergies are not the same, Eq. ( |3l| ) is not satisfied. In 
that case, Mn(N, H) should be calculated from Eq. (|27| ) 
with numerically solved jl(N,H) in Eq. (pi]). In order 
to confirm the above discussion, we show fl(N, H) and 
Mat (TV, H) in the cases of A Q = Ap and A a ^ A^ in 
the two-band system (J = a, [3). We set m a — 0.5mo, 
mp = mo, g a — gp = 0, T — 0.0001/io, A^ = and 
A a = 0, 0.34/^o and 0.68/iq. By assuming the parabolic 
band, pj — mjjl-K. The sum over r in Eqs. ( p4] ) and ( |27| ) 
is taken up to 30. We have checked that the results are 
the same within the numerical errors when we take the 
sum over r up to 50. It can be clearly seen from Figs. |l| 



that p,(N,H) and Mn(N, H) oscillate as the same func- 
tion of H when A a = Ap. When A Q is not equal to Ap, 
p,(N,H) is different from Mn{N,H), as seen in Figs. |2| 
and HI 
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FIG. 1. ft(N,H) (dotted line) and M N (N,H) (solid line) 
as a function of fxo/u)g when A a ~ Ap = 0. 
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FIG. 2. The same as Fig.[j] when A a = 0.34^ and = 0. 
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FIG. 3. The same as Fig.fi] when A a = 0.68/x and = 0. 
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In Fig. H we show M N (N,H), M^(p,H) and their 
Fourier transform intensities (FTIs), where we set m a — 
0.5m , trip = m , A Q = 0.68^ , A^ = 0, g a = gp = 
and T = 0.0001/io- With these parameters the fre- 
quencies are f a = m a (po — A a )/e = 0.16mo/xo/e and 
//3 = mo/io/ e - Thgse. parameters are the same as in 
the previous studiesErlia except for the temperature; the 
ground state energy is calculated directly by filling Lan- 
dau levels at T = in the previous studies. It is seen from 
Fig. H that the FTI has large peaks at the combination 
frequencies fp+f a , fp + 2/ a) 2/g + 2f a and fp - f a be- 
sides at the fundamental frequencies fp and /n^i-which are 
in agreement with previous results at T = OffiLl These 
combination oscillations are known as the quantum in- 
terference oscillations. The quantum interference oscilla- 
tions are caused by the oscillation of p. If p, is fixed, as 
seen in dotted line in Fig. U the FTI of M^(p,H) has 
peaks at the fundamental frequencies (fp and f a ) and 
higher harmonics (2fp and 2f a ) and it has no peaks at 
the combination frequencies. 
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given by R\ - r . Although the temperature dependence 
of the FTI at / Q is well fitted by - r , the temperature 
dependences of the FTIs at fp, fp± f a are not fitted 
by the temperature reduction factors, j r - Therefore, 
we should be careful to determine the cyclotron mass 
from the temperature dependence of the FTIs in two- 
dimensional systems. 
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FIG. 5. T-dependence of the FTIs of the magnetization at 
fa, f/3, fp — fa and fp + f a in the fixed N case (solid lines) 
and FTIs at f a and fp in the fixed fi case (dotted lines). 



B. special case: single-band system at T = 

In the single-band (pj — p, ujj = u), mj = m and 
Aj = A) and no-spin-splitting case (gj = 0) at T — 0, we 
can obtain Mj?(N, H) analytically. By using the identity 



(32) 



275 



[x] — x — — + — sin (2nrx) . 



where [x] is the largest integer satisfying [x] < x, Eq. (|2 
at T = is written as 



FIG. 4. M N (N,H) and M^(p,H) (lower figures) as a 
function of (Xa/up and their Fourier transform intensi- 
ties (upper figures), where the frequency has the unit of 
mopo/e. The Fourier transform is performed in the region 
of 255 < po/ up < 273. 

We show the T-dependence of the FTIs of M^{p,H) 



and M N (N,H) in Fig. 
same as those in Fig. 



, where the parameters are the 
In order to take in the sharp- 
ness of the oscillation at the ver y lo w temperature, we 
set the upper limit of r in Eqs. (^4|) and fl27|) as 75 be- 
low T/p = 0.00003. The T-dependence of the FTI is 
given by Rrj r (Eq- dH)) when p is fixed. When N is 
fixed, the temperature dependence of FTIs are no longer 



N(p, H) = 2 P uj 



p-A 1 
h - 

u 2 



(33) 



When we take this equation as a equation giving p(N, H), 
we get 



p(N,H) 



N 



2pu> 



1 



(34) 



This equation means that p is pinned at the Landau level 
and it jumps from one Landau level to the other when 
N/{2pLo) is integer. Since (p(N,H) - A)/lo + 1/2 is in- 
teger, we obtain from Eq. ( |l4"| ) that tl(p(N, H). H) does 
not oscillate as a function of H in the single-band case at 
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T = when N is fixed. By using the identity (Eq. (J3 
in Eq. we get 



u(N, H) = — + A + V — sin f 27rr ^ 

r— 1 x 



2p 



if 7 ' 



(35) 



where / = Nm/(2pe). 

It should be noticed that if Eq. ( p3] ) is considered as 
an equation giving fi(N,H) in the fixed N case, this 
equation is not satisfied when Eq. ( |34| ) is inserted: If 
we insert Eq. ([mJ) in Eq. (^3|), we get the right hand 
side as 2p(p(N, if) — A) for the single-band case, because 
{p(N,H) - A)/cj + 1/2 is integer in Eq. (§J). Although 
n(N, if )oscillates as a function of if, the left hand side 
of Eq. ( p3[ ) should be independent of if when N is fixed. 
Therefore, we cannot insert Eq. ( |34| ) obtained at T = 
into Eq. (p3|). The reason for this paradoxical result is 
the following. The Landau level with the energy fi(N, if) 
is partially filled at T — 0. This information is lost in 
Eq. ( jjji]) . In the same way we cannot get Mn [N, if) if 
we first set T = to get /Lt(iV, if) and insert it to get 
M N (N,H) = M^fi(N,H),H) in Eq. ©. 

The relation Eq. ([u]) is always satisfied in the single- 
band case even at T = 0. From Eqs. ( |3l| ) and j35|), we 
get 



Mn(N, H) = —Y—sm( 27rr4 
m z — ' nr \ H 

r=l v 



(36) 



We see from Eq. ( |2l| ) and Eq. (|36| ) that the oscillation 
of the magnetization for the fixed N has a phase-shifted 
and sign-changed saw-tooth shape comparing with that 
obtained in the fixed \i in the two-dimensional single- 
band spinless case at T = 0. 



IV. QUASI-TWO-DIMENSIONAL 
TIGHT-BINDING ELECTRONS 

In this section, we study the crossover from 
two-dimension to three-dimension in the quasi-two- 
dimensional electrons at T = 0. If we do not take into 
account the periodicity of the system, the system has 
the closed Fermi surface even for the very large effec- 
tive mass in the z-direction. Therefore, in order to study 
the crossover from two-dimension (with cylindrical Fermi 
surface) to three-dimension we have to study the periodic 
system. 

We take the tight-binding electrons in the simple or- 
thogonal lattice with lattice constants a, b and c and the 
hoppings t x , ty and t z . For simplicity we neglect the spin 
(i.e. g = 0) and set t x = t y = t. The Hamiltonian in the 



absence of the magnetic field is given by 



Wo = 5>t(k)£(k)c(k), 



(37) 



£(k) = — 2t(cos ak x + cos bk y ) — 2t z cos ck z 



(38) 



The magnetic field is introduced by the Pcicrcs substitu- 
tion, 



k + eA, 



(39) 



where A is a vector potential. In this paper the magnetic 
field is taken to be along the z-direction and the Landau 
gauge is used, 

A= (Hy, 0,0). (40) 
Then the Hamiltonian is written as 

H = — t {exp(iafc x )c^ (k x , k y — S, k z )c(k) + h.c.} 

k 

+ ^2 (-2tcos(bk y ) - 2t z cosfc z )c f (k)c(k). (41) 



where 



5 = 



eaH 



2n 



hc 4> Q b 



(42) 



4> = abH is the flux passing through a unit cell, and <f>o = 
2irhco /e is the unit flux quantum. For the numerical 
calculation we take the number of the flux quantum per 
plaquette, 



h = 



(43) 



to be a rational number. We use h instead of H here- 
after. If we take a = b w 10 [A] (the typical values 
for quasi- two-dimensional organic conductors), h 1/40 
corresponds to H = 100 [T]. We obtain the eigenvalues 
(ti(h)) by diagonalizing Ti. numerically. Then we calcu- 
late the Helmholtz free energy F(N, h) at T — 0, which 
is just the total energy, as 



1 N 

F(N,h) = w J2 e ^> 



(44) 



where N s is the site number. We fix the electron number 
to be 1/3 filled, i.e. N/N s = 1/3. The Fermi surface at 
h = and t z — is shown in Fig. |(| The area of the 
closed Fermi surface is 

/ = g/BZ, 

where /bz = 4,n 2 /ab. The magnetization is obtained 
from the numerical differential 



M N (N, h) 



5F{N, h) 
Sh ' 



where 



The chemical potential is given by 
fi(N, h) = e N (h). 



(45) 



(46) 
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FIG. 6. Fermi surface of the two-dimensional tight-binding 
electrons with filling factor 1/3. 

In the tight-binding model, even in the two- 
dimensional case (t z = 0), the density of states is not 
the sum of the delta functions but it consists of mini- 
bands with the finite width (Harper broadening). How- 
ever, the. width is very small especially in the case of 
t a = tf,E3 as shown in Fig. 0. On the other hand, the 
mini-gap, which means the spacing between the Landau 
levels, is of the order of 8th, as indicated by dotted lines 
in Fig. [?]. (The total band width is of the order of 8t, 
which consists of q delta-function-like mini-bands when 
h = p/q with co-prime integers p and q. However, p mini- 
bands are very close each other, making one super-mini- 
band. Therefore, the mini-gap is roughly of the order 
of 8t/(q/p) — 8th. For example, consider h — 1/40 and 
h = 3/118. These two cases should have similar density 
of states as can be seen in Fig. [t]. Indeed, they have 
the similar mini-band-gap, 8t x (1/40) « 8f x (3/118), 
if we neglect the mini-mini-gap inside the super-mini- 
band which consists of three mini-bands in the case of 
h = 3/118.) 

When t z becomes finite, the density of states becomes 
broader as shown in Fig. |^. The width of each mini-band 
becomes the order of At z . If the width of the mini-band 
is the same order of the mini-gap at t z = 0, that is, 
At z w 8th, the density of states has no energy-gap as 
in the three-dimensional systems in the magnetic field. 
We expect the crossover from two-dimension to three- 
dimension occurs at h ss t z /2t. In fact, we can see this 
crossover from the density of states for t z /t = 0.04 and 
l/h = 40, as shown in Fig. ||. The actual dimensional 
crossover in the quasi-two-dimensional organic conduc- 
tors may occur at H m 20 [T]. As, in these organic con- 
ductors, 8i.is.of the order of 1 [eV] and t z is of the order 
of 10 [K]B'E3, the crossover field is h rj t z /2t » 1/200, 
that is, ffw20[T]. 
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FIG. 7. Density of states for the tight-binding electrons in 
the magnetic field in the two-dimensional system (t z — 0). 
The arrows show the chemical potential for 1/3-filled band. 




FIG. 8. Density of states for the tight-binding electrons 
in the magnetic field in the quasi-two-dimensional system 
(tz/t = 0.04). The arrows show the chemical potential for 
1/3-filled band. 



7 



When the system is two-dimensional (t z = 0), the 
chemical potential is pinned in the narrow region of the 
non-zero density of states as in the case of the two- 
dimensional free electrons. In this case, the chemical 
potential and the magnetization oscillate periodically as 
a function of 1/h, as shown by the dotted lines in Fig. |j. 
The frequency of the oscillation of li(N, h) and Mn{N, h) 
is 1/3, which is in agreement with the closed Fermi sur- 
face area in the unit of the Brillouin zone, ///bz- 

In the quasi-two-dimensional system, the amplitudes of 
the oscillations of li(N, h) and Mjy(N, h) become smaller 
as t z increases as seen in Fig. ^. In particular, the os- 
cillation of n(N, h) is very small when the density of 
states becomes gapless ( h — 1/40 and t z /t = 0.04 as 
shown in Fig. |8J, where At z is nearly equal to 8th, i.e., 
4± z /(8th) = t z /(2th)=0.8). 



-0.3 



to 



-0.5 




2 4 3 2 yjj 4 



FIG. 9. jj, (upper figure) and Mn{N,K) (lower figure) as a 
function of 1/h. Arrows indicate the magnetic field for which 
the density of states are plotted in Figs. (?) and 

We compare the magnetizations in the fixed N and 
fixed \x cases in Fig. [l^. The magnetization with fixed ll 
at T = is obtained from the numerical differential of 
the thermodynamic potential Q(fi,h) as 



M^,h) 



sn(n,h) 

6h ' 



where 



n(Lx,h) = ^-J2^-^- 



(47) 



(48) 



The chemical potential is fixed to be u/2t = —0.414 to 



the magnetiza- 
30th cases of fixed 



satisfy N/N s = 1/3 at h = 0. In Fig. 
tion is plotted as a function of 1/h for 
\x and fixed N. We calculate Mn{N, h) and M^fj,, h) for 
two-dimensional case (t z = 0) and quasi-two-dimensional 



case (t z /t = 0.04). As expected, for the two-dimensional 
case (t z = 0), two curves have inverted and shifted saw- 
tooth shape each other, which can be seen from the upper 
figures in Fig. [To]. 

In the quasi- two-dimensional case, the difference be- 
tween Mn(N,K) and M M (/i, h) can be seen at h > 1/24 
and it becomes small at smaller field as clearly seen in 
Fig. [nj. This is consistent with the behavior of n(N, h) 
shown in Fig.^J, i.e. the amplitude of the oscillation of 
n(N, h) becomes small when t z / (2f) > h. 




FIG. 10. M N (N,h) and M M (ju, h) as a function of 1/h. 
Arrows indicate the magnetic field for which the density of 
states are plotted in Figs, fj] and ^. 



V. QUANTUM INTERFERENCE OSCILLATION 
IN QUASI-TWO-DIMENSIONAL ELECTRONS 

In this section we study the dHvA oscillation in 
the quasi-two-dimensional systems which has the small 
cylindrical Fermi surface and the quasi-one-dimcnsional 
Fermi surface, as shown in Fig. O. This Fermi sur- 
face is realized in k-(BEDT-TTF) 2 Cu(NCS) 2 , a-(BEDT- 
TTF) 2 KHg(SCN) 4 and GaAs/AlAs heterointerface. The 
small cylindrical Fermi surface is separated from the 
quasi-one-dimensional Fermi surface on the first Brillouin 
zone edge (dotted lines in Fig. [ll]) by the periodic poten- 
tial. In the semiclassical picture, electrons are expected 
to execute a large closed orbital motion by tunneling the 
first Brillouin zone gap at high fields. Then, the oscilla- 
tion due to the large closed orbit appears in the dHvA os- 
cillation. This is known as the magnetic breakdown phe- 
nomena. In k-(BEDT-TTF) 2 Cu(NCS) 2 and a-(BEDT- 
TTF) 2 KHg(SCN)4, the magnetic breakdown is observed. 
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FIG. 11. Fermi surface of the two-dimensional 
tight-binding electrons in the presence of periodic potential 
with filling factor 7/18. The area of the electron pocket cen- 
tered at (0, ±7r/2) is fee/ /bz ~ 1/18 in the unit of the area 
of the extended Brillouin zone ((2tt) 2 ). The orbit enabled by 
the magnetic breakdown has the area fp/fwi ~ 7/18. 

In order to obtain this Fermi surface we add the peri- 
odic potential 

V = V (-l) % c f (r)c(r) 

r— (ai x ,bi y ,ci z ) 

= y^{c t (k + g)c(k)+h.c.}, (49) 

k 

where 

g=(0,|,0), (50) 

to the Hamiltonian Eq. @ (when h = 0) or Eq. @ 
(when h + 0). When V ^ and h = 0, the Brillouin 
zone is folded, i.e. k x 6 [— 7r/o, 7r/a], fc y £ [— 7r/26, 7r/26] 
and fc^ G [— 7r/c, it/c]. In this section, we set electron 
filling factor N/N s to be 7/18 and V/2t = 0.1. In this 
case, fi/2t w —0.254 when /i = 0. In Fig. [ll] we plot the 
Fermi surface at t x = and /i = 0. The small electron 
pocket has the area 

fa ~ Tg^ BZ 

and the large orbit enabled by the magnetic breakdown 
at the gap on the zone edge has the area 



N fixed 
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FIG. 12. M N (N,h) and M^(n,h) (lower figure) as a func- 
tion of 1/h and their Fourier transform intensities (upper fig- 
ure) at tz — 0. The Fourier transform is performed in the 
region of 12 < 1/h < 62. The peaks of the FTI occur at f a , 
fp, fp — f a and fp + f a besides smaller peaks at 2f a , fp ± 2f a 
etc. 

We calculate the magnetization and its FTIs numer- 
ically as in the previous section. In Fig. [l2] we show 
the results in the two-dimensional case. It can be seen 
from the lower figure in Fig. [l^ that the amplitude of 
the oscillation with the frequency fp becomes large due 
to magnetic breakdown as 1/h becomes small. In the 
magnetic breakdown system, the oscillation with fp + f a 
is expected to exist but the oscillation with fp — f a is 
forbidden in the semiclassical picture. This semiclassi- 
cal picture is consistent with our quantum calculation as 
seen in the dotted liag in the upper figure in Fig. 
as known previouslyBQ. If the two bands are not con- 
nected with magnetic breakdown, combination frequen- 
cies (fp ± f a etc.) do not appear in M m (/j,, H), as can be 
seen from dotted lines in Fig. ^. On the other hand, in 
the two-band system, both frequencies fp ± f a are seen 
in M N (N,H) (Fig. |), where the FTI at fp + f a is larger 
than that at fp — f a - In the magnetic breakdown system, 
however, the FTI at fp + f a in M^(N,H) is very small 
as seen in the solid line in the upper figure ia-Eig. |lj. 
This behavior has also been known previouslyaQ. 

In the previous papers we have shown that the FTI at 
fp + f a in Mn (N, H) in the magnetic breakdown, pastern 
is anomalously enhanced by spin-splitting effecttJ'ta. We 
have interpreted the anomalous spin-enhancement of the 
FTI at fp + fa in Mjy(N, H) in the magnetic breakdown 
system as follows. From the results of two-band and mag- 
netic breakdown systems, it is found that the oscillation 
with fp + fa is caused by both effects of the magnetic 
breakdown and the chemical potential oscillation, as can 
be seen in the fixed fx case in Fig. O and in the fixed 
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N case in Fig. [|, respectively. If the electron spin is 
neglected in magnetic breakdown systems in the fixed N 
case, these two origins of the fp + f a oscillation almost 
cancel each other resulting in the small FTI at fp + f a 
in Mn(N,H). The spin-splitting may reduce the chem- 
ical potential oscillation. As a result, the cancellation 
becomes imperfect and the FTI at fp + f a in Mjv (N, H) 
in the magnetic breakdown system is enhanced due to 
the spin-splitting. 

The similar effect is expected in the case of quasi-two- 
dimensional system. As t z is increased, the chemical po- 
tential oscillation is suppressed, resulting in the enhance- 
ment of the FTI at fp + f a . This is indeed the case as we 
show in Figs. |lT 
and its FTIs for 



0.006 



and[14[ In Fig. |13| we plot M N (N,h) 
t z /t=0, 0.01, 0.02 and 0.03. As t z in- 
creases, the FTIs at f a , fp and fp — f a decrease as in 
the single-band case studied in the previous section. The 
FTI at fp + f a , however, depends very weakly on t z ; it 
gradually increases as t z increases, and it has a maximum 
around t z /t m 0.025 as shown by the arrow in Fig. [l4| . 




-0.002 



40 1/h 60 



FIG. 13. Mn(N, h) (lower figure) as a function of 1/h and 
their Fourier transform intensities (upper figure). The Fourier 
transform is performed in the region of 12 < 1/h < 62. The 
peaks of the FTI occur at f a , fp, fp — f a and fp + f a besides 
smaller peaks at 2f a , fp ± 2f a etc. 




0.01 tz/t - 02 



0.03 



FIG. 14. Main FTIs of M N (N,h) as a function of t z . An 
arrow indicates a maximum. 



VI. CONCLUSION 

We find that the dHvA oscillation in two-dimensional 
multi-band systems is given by Eq. (t27|) w ith the chemi- 
cal potential obtained by solving Eq. (|24|) . It is the gen- 
eral formula, i.e., it can be applied in the multi-band sys- 
tems at finite temperatures including the electron spin. 
In two-dimensional systems, the cyclotron effective mass 
is generally not obtained by fitting the temperature de- 
pendence with the reduction factor Rtjv It should be 



obtained from the fitting by Eq. (|27p with Eq. (|2J). 

We also study the crossover of the dHvA oscillation 
from two-dimensional systems to three-dimensional sys- 
tems. In the simple case of the single Fermi surface we 
show that the oscillation of the magnetization calculated 
in the canonical ensemble (fixed N) approaches smoothly 
to that calculated in the grand canonical ensemble (fixed 
/i). In the magnetic breakdown system we find that the 
Fourier transform intensity at fp + f a gradually increases 
as a three-dimensionality increases and it has a maxi- 
mum at t z /t w 0.025. The t z dependence can be ob- 
served experimentally by applying the uniaxial stress in 
the quasi-two-dimensional organic conductorsaCJ, since 
the uniaxial stress along the z-axis will increase t z . 
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